Coulomb Gas and Sine-Gordon Model in Arbitrary Dimension 
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Exact phase structure of the neutral Coulomb gas has been mapped out by the functional renor- 
malization group study of the equivalent sine-Gordon field theory in arbitrary dimension. The high 
temperature fixed point is obtained and it is shown that the Coulomb gas has a single phase for 
d > 2. Renormalization group equations of the 3-dimensional sine-Gordon model are compared to 
those of the vortex-loop gas i.e. the gas of topological defects of the 3-dimensional isotropic XY 
spin model which demonstrates that the two models belong to different universality classes. 
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I. INTRODUCTION 



Universality which expresses the idea that different mi- 
croscopic physics can give rise to the same scaling be- 
havior at a phase transition, is one of the most powerful 
tool of critical phenomena. For example, if two different 
models can be mapped onto each other then they should 
belong to the same universality class hence by the study 
of one model one can obtain the critical behavior of the 
other one. If the critical behaviors of two models agree 
or disagree it is possible to decide whether they belong to 
the same universality class. Furthermore, models of the 
same class of universality can be used to test and com- 
pare various types of methods suitable for the study of 
critical behavior. A canonical example for model equiva- 
lences is the sine-Gordon (SG) scalar field theory defined 
by the Euclidean action [l|-|4| 



S 



d d x 



-(d ll ip x ) 2 +ucos(/3(p x ) 



(1) 



which belongs to the^ universality class of the neutral 
Coulomb gas (CG) [j| and the isotropic XY classical 
spin model [1] for d = 2. They under go a Kosterlitz- 
Thouless-Berezinski [6[ (KTB) phase transition which 
was well reproduced in the framework of the O(N) sym- 
metric scalar field theory for d = 2 and iV = 2 
Furthermore, the SG scalar theory is the bose form of 
the 2-dimensional fermionic Thirring model Q and the 
bosonic counterpart of multi-flavor Quantum Electrody- 
namics and multi-color Quantum Chromodynamics are 
also SG type models Q. Although bosonization is well 
established only in d — 2 dimensions, it is known that 
the mapping between the SG theory jO} and the neutral 
CG holds in arbitrary dimension [l|, y| . 

The goal of this work is (i) to consider possible equiv- 
alences between CG, SG, XY models in d ^ 2 dimen- 
sions, (ii) to map out the exact phase structure of the 
d-dimensional neutral CG by the functional renormaliza- 
tion group (RG) study of the equivalent SG scalar theory. 



II. RG STUDY OF THE SG MODEL 

By the inclusion of the wave function renormalization 
in the functional RG study of the SG scalar field theory 
performed in arbitrary dimension we go beyond the ap- 
proximations used in previous works [2|-|j]. We use the 
effective average action functional RG method .10] where 
the evolution equation reads as 



kd k T k 



Tr[(R k + dfrTkltp])- 1 kd k R k ], (2) 



where k is the running momentum cutoff, the trace Tr 
stands for the integration over all momenta. We intro- 
duce the power-law regulator function R k = p 2 (fc 2 /p 2 ) 
where b > 1. Let us note, that the field carries a di- 
mension for d ^ 2 thus the frequency of the SG model 
([T]) becomes a dimensionful parameter for d ^ 2, i.e. 
2 = k 2 ~ d /3 2 where f3 is dimensionless. Eq. ([2} has 
been solved over the functional subspace spanned by the 
ansatz for the SG model 



r, = / d d 



i i 



y^{d^fx) 2 +Uk cos(fc~ 



fx) 



(3) 



where the local potential contains a single Fourier mode 
and the field has been rescaled as ip — > <p/f3. In case of 
the d-dimensional SG model it is convenient to introduce 
a dimensionless variable (p = kS 2 ~ d ^ 2 tp, then ([3]) reads as 



d d x 



Zk{d^ x ) + V k (<p x ) 



V k {fx) = u k cos((f x ) 



(4) 



where z k = 1/0$, = l/(fc 2 - d /3f) stands for the field- 
independent wave-function renormalization which has a 
dimension of k d ~ 2 . The RG equation has been modified 
if a dimensionless field variable ip has been introduced, 



kd k T k [tp] = -Tr 



kd k R k 



Rk + k d - 



2 nm 



(5) 



with the notation ' = d/dip. Let us note the term — [(d — 
2)/2](pd<f, is neglected on l.h.s. of (f5]) because it does not 
give contribution to the periodic part of the action Q 



2 



and the blocked potential V k ((p x ) contains only periodic 
interactions (and z k is field-independent). Although RG 
transformations generate higher harmonics, we use the 
ansatz (Q} which contains a single Fourier mode since 
in case of the 2-dimensional SG mod el it was found 
to be an appropriate approximation 12]. Then Eq. ([5]) 
leads to the evolution equations for the blocked potential 
Vk(<p) and the wave-function renormalization Zk, 



d k V k 
kd k z k 



d d p 
(2tt)° 

2tt 



V k kd k R k , 



2tt 



- Hip Uk 



(I 



'- 2 - d vl"? 



d d p 

(2n) d 
dp 2 



(6) 



(7) 



withP fc = (z k k 2 ~ d p 2 +R k +k 2 - d V k y 1 . Let us note since 
the l.h.s of ([7]) is independent of the field, a projection 
onto the field- independent subspace has been introduced 
on the r.h.s of 0. The scale k covers the momentum 
interval from the high-energy/ultraviolet (UV) cutoff A 
to zero. Inserting the ansatz (HJ) into Eqs. ((6|), (0, flow 
equations for the dimensionful couplings are 



kd k u k = 
kd k z k = 



kd k R k P- y/P 2 - (k 2 - d u k ) 2 



„2-d 



p 

kd k R k 



y/P 2 - (k 2 - d U k ) 2 



(8) 



-{k 2 - d U k ) 2 P{d p 2P + 



[p2 _ ( fc 2-d Ufc )2j5/2 



(k 



2-d 



u k ) 2 p 2 (d p 2P) 2 (4P 2 + (k 2 - d u k ) 2 ) 



,(9) 



where P 



d[P 2 - {k 2 - d u k ) 2 ]V 2 
2 ^R k and J = Jdpp d - 1 n d /(2n) d 



z k k- "p 

with the d-dimensional solid angle f^. Momentum in- 
tegrals have to be performed numerically, except the 
linearized form of Eqs. ©, © around the Gaussian 
fixed point where analytical results available. The lin- 
earized RG with dimensionless couplings u k 
k 2 ~ d z k reads as for d = 1 dimensions 



(1 + kd k )u k 
(-1 + kd k )z k 



1 



46sin(^) 



1 — 2b 

z k 2b Uk 



0(ul) 



Cl{b) „ 



0{u 



(10) 
(11) 



Finally, for d — 3, linearized RG equations are 

(3 + ^ fe K = ^^^^ + 0(u 2 ) (14) 

-^-z^ul + Oiul) (15) 



(1 + kd k )~z k 



with c 3 (6) = xg^^nj'it/t-jb)] ■ Due t0 the tree-level scaling 
of Zk, the non-trivial fixed point appears for d < 2 and 
disappears for d > 2 in the RG flow. However, a "turning 
point" can be identified for d > 2 where the irrelevant 
coupling u k turns to a relevant one. For d = 3 the turning 
point is at z* = [247r&sin(7r/(26)]^ 2b /( 1+2fc ). 

In general, RG equations ([5]) and © (exact for a single 
Fourier mode) have to be solved numerically. Important 
feature of the exact RG flow the emergence of a new low- 
energy /infrared (IR) fixed point related to the degener- 
acy of the blocked action. Namely, Eqs. jSJ, ((9]) become 
singular at the momentum scale where k — k 2 ~ d u k = 
with k = min p 2 P = bk 2 [S k /(b - l)] 1 ^ 1 / b . Therefore, it 
is convenient to redefine the dimensionless coupling con- 



stant as u k 



= u2-d 



u k /k = k u/k which tends to one in 



case of degeneracy. Exact RG equations ©, © were 
solved for b = 2 and flow diagrams are plotted in Fig.Q] 
for d = 1, and in Fig. [2] for d = 3 dimensions. The phase 




FIG. 1: Phase diagram of the SG and the equivalent CG 
models for d — 1 dimensions (similar RG flow can be drawn 
for d < 2) . Arrows indicate the direction of the flow. Param- 
eters u, z are related to the fugacity y and the temperature 
T. 



with 0,(6) = »^i^ffl 



Eqs. CEU]), IIITJ have a non- 
trivial fixed point at z ± = [46sin(7r/(26)] 2fc /( 1 - 2b ), = 
(4?r / ! ci)zt~ b ^ 2h \ For d = 2, linearized RG equations are 



(2 + kd k )u k = —-u k + 0{ut) 
Airz k 

C2(5) z^u 2 4 



kd k z k 



o{ul) 



(12) 
(13) 



with c 2 (b) = ^i^^l/b] an d result in a KTB type (i.e. 



infinite order) phase transition with z+ = l/(87r) [lllll 



structure is independent of b but we use b = 2 here since 
it was found to be an optimal renormalization scheme for 
SG type models [Tj3]. The IR fixed point corresponds to 
the degeneracy was found in any dimensions. However, 
the d = 1 case possibly needs further improvement since 
its phase structure shows similarities to the QM anhar- 
monic oscillator where spontaneous symmetry breaking 
is not allowed [l4j]. There the higher order terms in the 
gradient expansion handle the problem. 

Let us emphasize that only the exact RG flow with the 
inclusion of the wave function renormalization is suit- 
able for the determination of the IR fixed point since the 
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FIG. 2: Phase diagram of the SG and the equivalent CG 
models for d — 3 dimensions which indicates a single phase 
for d > 2. 



perturbative (truncated) RG equations are non-singular 
[2|. Thus any previous attempts based on perturbation 
theory [21 or by the usage of the local potential approx- 
imation J, Q were unable to determine the IR behav- 
ior in a reliable manner. Furthermore, the exact RG 
flow changes the position of the non-trivial fixed point 
obtained for d — 1. Indeed, for b = 2, linearized RG 
equations (JTUJ*, CP result in 1/z* = 2 10 / 3 « 10.08 and 
= u+/(2\fzl) = 2/vl5 w 0.52 which have been mod- 
ified by the exact RG, see Fig.[TJ Similarly, the turn- 
ing point of the d = 3 case determined by the exact 
RG coincides with that of the linearized RG which is 
1/z* = 4(3tt) 2 / 3 « 41.9 only for small u, see Fig.[2J Fi- 
nally, let us mention that the comparison of exact and 
perturbative (linearized) RG flow can be used for opti- 
mizing the renormalization schemes. For d = 1 and 6=1 
the non-trivial fixed point determined by the linearized 
RG (1/z* = 16 and = = 4/v / 5 ~ 1.79) is out of the 
range of validity of the exact flow (Uk — Wfe < 1), hence 
it represents a drawback on the applicability of the RG 
scheme b — 1 for d ^ 2 dimensions. 



III. PHASE STRUCTURE OF THE NEUTRAL 

CG 



case of low densities (limit of low fugacity) and there are 
different competing views of what happens to the neutral 
plasma at high densities [H, [17|. Thus, an RG study of 
the d-dimensional CG beyond the dilute gas approxima- 
tion is required which is the goal of the present work. 

Since the mapping between the CG and SG models 
holds in arbitrary dimension [3J, |5[ (and it is exact in case 
of point-like charges) the RG study of the d-dimensional 
SG model can be directly used to map out the phase 
structure of the CG model. In the framework of pertur- 
bation theory (limit of low fugacity) RG equations of the 
CG model in arbitrary dimension [l5| reads as 



dx 



-x(xy 2 +d- 2), 



dy 
dl 



-y{x-d) (16) 



with di = —kdt, y ~ Uk and x ~ 1/T ~ l/zf. where 
T is the temperature and y is the fugacity. Using these 
identities Eq. (Ti"6l) can be rewritten as 

1 



((d-2) + kd k )z k 



-c z u k , (d+kd k )u k = c u —u k , (17) 

Zk 



with constants c u , c z and they are found to be similar to 
the linearized RG equations obtained for the SG model, 
e.g. the non-trivial fixed point of the 1-dimensional SG 
model can be identified in the flow generated by (|17l) . In 
the sharp cutoff limit (b — > oo) the linearized RG equa- 
tion derived for the Fourier amplitude of the SG model 
reduces to (d + kdk)u k — QdUk/(Sk2(2ir) d ) which is iden- 
tical to the second equation of (fTTj) . However, the sharp 
cutoff docs not support to the derivative expansion at 
higher order, i.e. for b — > oo the multiplicative constant 
c c i{b) of the RG equation obtained for z k becomes infinity. 

Most important new results of the exact RG flow are 
the existence of the high temperature (u* = 1, 1/5* = 0) 
and the absence of new further non-trivial fixed points. 
Let us note that, the perturbative (truncated) RG has 
a non-singular structure hence it is not able to recover 
the degenerate potential therefore it is not suitable to 
observe the high temperature fixed point [L5| . Moreover, 
the exact RG changes the position of the non-trivial fixed 
point for d < 2 and the turning point for d > 2. Further- 
more, since the mapping between the SG and CG models 
is exact for point-like charges the exact RG flow indicates 
a single phase for the CG for d > 2. 



The CG model has been the subject of an intense study 
in last decades @, H, H, ll"5T - fl"7| and there is a continuos 
interest in the usage of the SG representation of CG sys- 
tems @, H, [!, Oj} since it has an important relevance 
in various models of critical phenomena (e.g. superfluid 
transition in He films, melting of crystals in d = 2, arrays 
of Josephson junctions, roughening transition) . However, 
the complete understanding of the critical behavior of 
the d-dimensional neutral CG which is a plasma of equal 
number of positive and negative charges is still lacking. 
For example, previous investigations (even for the well 
studied 2-dimensional CG model) were restricted to the 



IV. ISOTROPIC CLASSICAL XY SPIN MODEL 

The partition function of the d-dimensional isotropic 
classical XY spin model |l| can be expressed in terms of 
topological excitations of the original degrees of freedom. 
For d = 2 dimensions, the dual theory is the vortex gas 
which is known to belong to the class of universality of 

the neutral CG [a, nam. Two-dimensional generalized 

models are well known where both the CG and the vortex 
gas are included as particular limiting cases [H, [3 H3 
and are self-dual under the duality transformation. For 



4 



d = 3. the dual theory is the gas of interacting vortex 
loops (IS H3 (i-e. the lattice CG HH). Corresponding 
flow equations have been derived for the parameters K 
(i.e. the coupling between the spins) and y (i.e. the 
fugacity of the vortex loops) by real-space RG method 
(a ~ 1/fc) in the limit of low fugacity [a, Sill, 

4.7T 3 

ad a K = K - —K 2 y, ad a y = (6 - tt 2 KL) y, (18) 

where L approaches a constant in the IR limit (L — > 1) 
20], or it is weakly divergent (L —> ln(a/a c ) + 1) [2lj . 
Since ad a — —kdf, and by using the identities y = u, 
K = 1/5, RG equations of (jT5J) are rewritten as 

47T 3 / 1 ~K 2 L \ 

kd k h = h 3~"fe' kdkUk = ( Y~2 3 ) " fe ' 

which can be compared to the linearized RG obtained 
for the SG (P | . (fT5 ]l and the CG dTTJ) models in d = 3 
dimensions. The only disagreement is the sign of the tree- 
level scaling term of Sk but it has important consequences 
since the RG flow of the vortex-loop gas has a non-trivial 
fixed point which is absent in case of the 3-dimensional 
SG and CG models. Therefore, it demonstrates that the 
vortex-loop gas has a different critical behavior, thus it 
belongs to a class of universality differs from that of the 
SG and CG models for d ^ 2 (for d = 2, couplings K, z 
have no tree-level scaling thus SG, CG, XY models are in 
the same universality class). Let us emphasize that RG 
equations are compared at the linearized level, however 
the exact RG study of the SG and CG models is required 
since it shows the absence of new non-trivial fixed points, 
thus there is no room to find a mapping between the 
parameters of the vortex-loop gas and the 3-dimensional 
CG which could produce the same phase structure. 



V. SUMMARY 

In this work, the exact phase structure of the neutral 
CG was mapped out by the functional RG study of the 
equivalent SG scalar field theory in arbitrary dimension. 
By the inclusion of the wave function renormalization 
and the determination of the exact RG flow for a sin- 
gle Fourier mode, we went beyond the approximations 
used previously. The high temperature fixed point of the 
d-dimensional CG was determined and shown that the 
perturbative RG is not suitable to recover it. The po- 
sition of the non-trivial fixed point of the RG flow for 
d < 2 and the turning point for d > 2 dimensions given 
by the exact flow is compared to the linearized one. The 
non-trivial fixed point of the CG for d < 2 shows similar- 
ities to the Wilson-Fisher fixed point of the O(N) model. 
The absence of new further non-trivial fixed points were 
also demonstrated which results in a single phase for the 
neutral CG with point- like charges for d > 2. This has 
two consequences, (i) no further phase transition can be 
identified at high densities for the neutral CG as opposed 
to e.g. [lj| [lTf, (ii) there is no way to find a mapping 
between the phase structure of the neutral 3-dimensional 
CG and the vortex-loop gas which is the gas of topologi- 
cal defects of the isotropic XY spin model for d = 3, thus, 
it was demonstrated they belong to different universality 
classes. 
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